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Abstract: The nonuniform black strings branch, which emerges from the critical Gregory- 
Laflamme string, is numerically constructed in dimensions 6 < D < 11 and extended into 
the strongly non-linear regime. All the solutions are more massive and less entropic than the 
marginal string. We find the asymptotic values of the mass, the entropy and other physical 
variables in the limit of large horizon deformations. By explicit metric comparison we verify 
that the local geometry around the "waist" of our most nonuniform solutions is cone-like 
with less than 10% deviation. We find evidence that in this regime the characteristic length 
scale has a power-law dependence on a parameter along the branch of the solutions, and 
estimate the critical exponent. 
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1. Introduction and summary 

The surprising discovery of Gregory and Laflamme (GL) ||l| that below certain mass the 
uniform black strings, existing in the backgrounds with compact extra dimensions, are 
perturbatively unstable created considerable excitement. It initiated intensive research 
aimed to understand the consequences of this instability and to figure out its endpoint. 
In this regard, some progress has been made in the numerical simulations of the decaying 
string Yet, the endpoint remains elusive as the simulations crash (due to numerical 
problems) when the spacetime is still highly dynamical and far from settling down to the 
static solution. A different program ||3| to determine the endstate proposes to find a phase 
diagram accommodating all possible static solutions. The hope is that in principle the 
endstate can be inferred from the diagram since by definition it is a static solution. This 
research has recently culminated in the reviews [Q, 

A proposed phase diagram is portrayed in figure |l|. The vertical axis shows the di- 
mensionless mass density while the horizontal axis is related to the scalar charge. The 
branches of solutions represented by solid lines are presumably stable and the dashed lines 
indicate unstable solutions. The diagram has several remarkable features: 1) The merger 
point where the black hole branch meets the black string branch and where the horizon 
topology changes. 2) A critical dimension order = 13 above which the order of the 
phase transition changes from the first to the second Q and where the emergent branch of 
nonuniform solutions becomes (thermodynamically) stable Q. According to the diagram. 
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Figure 1: A suggested phase diagram. The vertical axis is the dimensionless mass density and 
the horizontal axis is related to the scalar charge. Shown is the GL point where the uniform string 
becomes marginally unstable and from which a new branch of non-uniform strings emerges. This 
branch extends until it meets the caged black hole branch at the "merger point" . In D > 13 
the NUBSs branch is less massive than the critical string indicating that the order of the phase 
transition, triggered by appearance of the GL-mode, changes from first to second. 



the endpoint of the GL instability depends on the dimension and it is either a caged black 
hole in D < 13 or a nonuniform black string otherwise. For now, only certain parts of the 
diagram in certain dimensions were verified by analytical/numerical constructions. This is 
described next. 

The black hole branch has been constructed numerically in five ^ and six ||9|, p!o| 
dimensions. In particular, it has been shown to approach the merger point, where the north 
and the south poles of the caged black hole are expected to intersect across the compact 
dimension [^. In other dimensions only the beginning of the BHs branch was constructed 
perturbatively from the Schwarzschild-Tangherlini [11| solution in [ll2| , [l^ , |l^ , [T5|| and 
numerically in ||l^. The branch of the NUBS was found perturbatively from the critical 



GL string in five |17], six |15] and in other, up to sixteen, dimensions [^]. In six |15] and 
in five [19| dimensions the branch was numerically extended into non-linear regime, where 
it was shown to behave in accordance with diagram |^. 

In this paper we aim to construct the NUBS branch in higher dimensions and to extend 
it beyond linear perturbations into the strongly non-uniform regime. To define the "non- 
linear" , "strongly non-uniform" regime we use a simple geometrical measure of the horizon 
deformation p^ , 

A-^ft^-l), (1.1) 

where Rmin/Rmax are the minimal/maximal areal radii of z = const sections of the horizon. 
A = for uniform strings and A ^ co at the pinch-off.^ It turns out that while at pinch- 



^ Apparently, Rmax remains finite in this limit. 
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off, Rmin 0, many other thermodynamical and geometrical quantities, describing the 
NUBSs, asymptote to finite values. This fact suggests a natural definition of the strongly 
non-linear regime as a situation when these variables are close to saturation.^ 

Recently, intriguing proposals regarding the local geometry near the "waist" of the 
extremely non-uniform strings have been put forward by Kol |^, 20 1. Specifically, in Q 



Kol advocates that locally the merger spacetime is cone-like, and moreover in I? < 10 the 
local slightly off-merger metric and its functions have power-law dependence on parametric 
distance from the critical (merger) cone. In other words, if p is a parametrization of the 
static NUBSs branch, such that p* corresponds to the merger, than quantities with the 
length dimensions are conjectured to scale as 

log(^) = 7 log(p — p*) + (p[y logijp — p*) + const] + const, (1-2) 

where (/> is a periodic function.'^ Namely the leading behavior is a power-law with an 
exponent 7; and this is dressed with periodic "wiggles". 

Unfortunately, the data available in literature provides only limited insights on what 
happens at the merger since neither along NUBS nor along the BH branches is the merger 
point approached close enough. Still, some evidence in favor of the local cone geometry 



was reported in |21|. The scaling (|l.2| ) has not been tested so far. 

In this paper, we partially close the gaps. We use numerical techniques to find fully 
non-linear branches of NUBSs in dimensions from six to eleven and confirm that these are 
in accordance with diagram ||. In dimensions 6, 7, 8, and 9 we are able to extend the 
NUBSs branch into the deeply non-linear regime. In this range of the dimensions we find 
the limiting values of the thermodynamic and geometric variables, these are listed in table 
|l[ Furthermore, since our method allows approaching the merger very closely, we compare 
the local near "waist" geometry with a cone. This is explicitly done in 6D where we have 
our most nonuniform solutions with the highest resolution. We find that the metrics agree 
well and the agreement improves with growing A. For our most nonuniform solution the 
discrepancy between the near-waist metric and the cone does not exceed 10%. 

Next, we verify the scaling (|1.2| ) and find evidence for the power-law. We extract the 
approximate critical exponents and list them in table ^. However, the wiggles are not seen 
for our non-uniform solutions. We believe that in order to actually test for the presence of 
wiggles, one must have solutions spanning over several orders of magnitude in A. Currently, 
this challenging task stands beyond the capabilities of our numerical experiments. 

In the next section we describe the setup and derive the field equations and the bound- 
ary conditions. We define the physical variables and describe how the cone and the scaling 
( |1.2| ) are tested. In section |3| we elaborate on our numerical method. In section ^ we 
describe the solutions and their properties and discuss them in section Various technical 
details are found in the appendices. 



^The saturation is best noted using A, since in A the measurables approach their limiting values roughly 
exponentially. 

^This sort of behavior is known to appear p2| p3| in the near-critical collapse of a scalar field (see 
| p5| ] for a D-dimensional version). Indeed, there are speculations that the merger cone might be related 
to the critical collapse by a double analytic continuation and change of the boundary conditions. 
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2. Static black strings 



We consider D-dimensional background of the form M^~^'^ x S^, that has one spatial 
dimension, z, compactified on a circle of length L, z ^ z + L. We are interested in static 
black objects, which are spherically symmetric in the extended dimensions. The simplest 
such a solution is the uniform black string whose metric in Schwarzschild coordinates is 
given by 

dsl = dsg^i^^^ + dz'^, 

^4cw, = -fip)dt' + fip)-'dp^ + p^dnj_„ (2.1) 

where f{p) = 1 — {po/p)'^~^, po designates the horizon location and dQ,'^_2 is the metric 
on a unit sphere §'^~^. Gregory and Laflamme showed that this string is unstable below 
certain mass [Q]. Here, we construct the branch of static non-uniform black strings that 
emanates from the GL-point. 

The most general metric describing the black strings can be written as 

ds2 ^ -e'^^dt^ + da{r, zf + e'^^'dn^-s, (2-2) 

where A and C are functions of r and z, da{r, z)^ is a two-dimensional metric in the (r, z) 
plane. The horizon is located where exp(yl) = 0. A convenient ansatz that we employ in 
this paper^ uses "conformal" coordinates in which 

ds^ = -f{p(r))e^^ + e^^'idr^ + dz^) + p{rfe^^ d^l_^. (2.3) 

When A = B = C = Q we reproduce the UBS solution in conformal coordinates (the 
transformation p ^ r from Schwarzschild to the conformal coordinates is described in 
appendix ^.) The horizon is located at the radius tq, where /(rg) = 0. Without loss of 
generality we assume rg = 0. 

2.1 Equations and boundary conditions 

In order to find the black strings we solve the Einstein equations, Rab = 0, that in our 
coordinates ( ^!^ ) split into three elliptic equations 

A ^ + [drAf + {d,Af + (£> - 3) {drAdrC + d,Ad,C) + 

+ ^^^^{2drA + {D -'i)drC) + ^ {2drA-{D -'i){D -A)drC) = 0, 
2py/f 2p 

AB - (2 drA + {D- ^)drC) dr-C - (2 d,A + {D - 4)0,0) d,C - 

- %rc - (2 oIa + (D - mo - 



2pVJ 2p 
(I)-3)(D-4) l-e2^-2^ 



0, (2.4) 



2 p2 

AC + {drA + {D- 3)drC) drC + {d,A + {D - 3)5^C) d,C + 

+ ^:^^drC + ^ {2drA + {3D - 8)drC) + {D-4) ^ = 0' 

2 py/f 2p p^ 



*We found that the expUcit separation of the UBS background is essential for stabihty of our numerical 
scheme. 
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where A = 9^ + 9^, and two hyperbolic constraints 



Gl : drzA + {D- 'i)drzC - dr-A {d,B - d,A) - 

- {D- ?,)drC {d,B - d^C) - drB {d^A + {D - ?,)d,C) + 



dlA - dlA + {D-3) [plC - die) + {d^Af - (drAf + 



+ (Z)-3) (d,Cf-{drCf 



2 d^Ad^B + 2dr AdrB 



2{D-3) id,Bd,C - drBdrC) 
Vf 



D-A 

777 



{drA - drB) + 



+ —[(D- 4)drA + (D- 2)drB -2(D- 3)drC] 
P 



0. 



(2.5) 



As first observed by Wiseman |26|, the Bianchi identities, G^.^ = 0, imply that the con- 
straints satisfy Cauchy-Riemann relations, 



dM + drV = 0, 
U = V^G:., V = 



drU - d,V 
1 



0, 



9{g:-gi). 



(2.6) 



This in turn means that each one of U and V separately satisfy the Laplace equation. 
Hence, the "constraint rule" follows: provided one of the constraints vanishes along the 
boundary of the domain and the second one is zero at a single point of the boundary — both 
constraints are guaranteed to vanish inside the domain [26|. 

The elliptic equations ( |2.4D are solved outside of the horizon. The periodicity in z and 
the reflection^ symmetry at z = suggest the sufficient domain of integration {{r, z) : < 
z < L/2,r > 0}. 

The equations are subject to boundary conditions (b.c.) The reflection symmetry at 
z = and the periodicity at z = L/2 are translated into 



z=0 



'z^\z=L/2 



for ^' = A,B,C. 



A regularity of (j^^ at the horizon supplies two conditions 

drA\r=0 = drC\r=0 = 0. 



(2.7) 



(2.8) 



However, we need three b.c. for our three fields. The missing condition is obtained from the 
"constraint rule", following which we impose U = along all boundaries. The constraint is 
automatically satisfied along z = and z = L/2 due to (2.7), and it vanishes exponentially 
at large r (as does any z-dependence in a Kaluza-Klein background). A regularity of U at 
the horizon yields 

(dzB - dzA) \r=o = 0, (2.9) 



^An appearance of the GL zero-mode spontaneously breaks the translational invariance along z. The 
remaining symmetry to the translation by half a period implies that the branch of NUBS solutions that 
emerges from the GL-point has reflection symmetry about z — 0. 
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and implies that the surface gravity is constant along static horizon. Integrating (|2.9| ) from 
z = L/2 we obtain a Direchlet boundary condition for i?, 



5(0, Z) = BU=L/2 + ^(0, Z) - 



(2.10) 



The integration constant Bq = -B| 2=^/2 is freely specifiable; once it is chosen the NUBS 
solution is unique for a given L. Bq = corresponds to a unform string, a positive Bq 
yields certain non-uniform string solution.^ 

Asymptotically we require that the spacetime is cylindrically-flat, A = B = C = 0. 



The leading fall-off of the fields in this region is |27| 

b 



A 



C ~ 



B ~ 



clog(r) 



for D 



C 



for L> > 5. 



(2.11) 



This form is useful in numerical implementations whenever the asymptotic boundary is 
placed at a finite r. 

Following the constraint rule and imposing V = at the asymptotic boundary will 
relate the asymptotic constants^ in ( ^.11 ). However, in our numerical method which has 
the outer boundary at a finite r^, we could get better overall accuracy by simply enforcing 
the fall-offs ( p. 11 ) along Tq. The constraints are verified in appendix^. 



2.2 Charges and geometry 

Thermodynamic properties of the solutions are determined by the asymptotic charges - 
the mass and the tension, and by the horizon measurables - the surface gravity and the 
horizon area. The variables are related by the generalized Smarr's formula. The geometry 
of the black-strings is conveniently visualized by embedding their horizon into flat space 
and by examining the proper length of the compact circle along r = const slices. 



Asymptotics and charges 

The charges are related to the asymptotic constants (|2.11|) by pi 



D-3 -1 
1 -(I) -3) 



m 




tL 





b 



(2.12) 



where i^o-i = D tt^^"^ /T{D/2 + 1) is the surface area of a unit sphere, and Pq~^/2 

originates from the background f{p{r)), see (|2.3D . (We set below po = 1). For convenience 
we use the "asymptotic adapted" units in which the d-dimensional Newton's constant 
Gd = Gd/L = 1. 



^An alternative way to specify the integration constant is to fix the surface gravity k. This determines 
B along the horizon as B — A — log k. We checked that both approaches produce comparable results. In 
either case, there is only one constant (either Bq or k) that has to be chosen. We fix Bo and in this case 
the thermodynamical k becomes a derived quantity. 

^To this end, one must also include the subleading term in the expansion of C, that for D > 5 is c/r^. 
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The dimensionless charges, defined with respect to the critical uniform black string's 
values (which can be read from (2.12) after setting a = b = 0), are 



rh = m/rric = 1 + 2 a - 2 b/{D - 3), 

f EE t/tc = 1 + 2a - 26(D - 3). (2.13) 

Notice the factor D — 3 that in higher dimensions increases the significance of b in the 
tension computation and suppresses it relatively to a in the mass calculation. 

One can also define the relative tension, n = rL/m. Normalizing with respect to 
nc = l/(-D — 3) one has 

l + 2a-2b(D -3) , 
n = n Uc = 7 n ^ t (2-14) 



Horizon variables 

The normalized temperature (or surface gravity) of the horizon reads 

f = k = k/kc = e^-^|^=o, (2.15) 

where T = k/ {2-k) and Kc = ^f'{po) = {D — 4)/2. The normalized horizon area is given by 



^ Jo 



(2.16) 



where Ac = LQ £)_3 is the surface area of the critical string. 

The horizon measurables are related to the asymptotic charges by the generalized 



Smarr's formula |27, 



{D-3)m = {D-2)ST + TL, (2.17) 

where S = A/AGd is the Bekenstein-Hawking entropy. 

The first law of thermodynamics relates variations of the thermodynamic variables. 
If the asymptotic length of the compact circle is kept fixed (as in our case), the first law 
reduces to 

dm = TdS dm = TdS dm = TdS . (2.18) 

rric D-3 ^ ^ 

Geometry 

We examine the black string's geometry by embedding the r = const hypersurfaces 
into flat space, 

dsji,, = dz^^ + Rldnl^,, 

dZe = e^("'") dz, 

= e^('^'^)p(r). (2.19) 
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The dimensionless proper length of the compact circle along an r = const hypersurface is 
given by 



Lpropir) dZ,{r) = i ^ e^^'-'^) dz. (2.20) 

In particular, the dimensionless proper length of the horizon is Lhor = ip-rop(O). The 
horizon's areal radius, Re{z) = e^'^'^'^\ is not constant along the z-direction. 

2.3 Near the merger point 

For progressively non-uniform strings the minimal areal radius of the horizon, Rmim is 
steadily shrinking. Kol argues that in the regime when Rmin — > the local geometry 
in vicinity of the "waist" becomes cone-like (see also additional arguments in recent [pg|].) 
Moreover, there is the power- low scaling ( |1.2[ ), which in D < 10 is dressed with periodic 
"wiggles" 0. 



A direct comparison of the 6D solutions of |18] with the cone metric was already carried 
in |2^, where in spite of low numerical resolution authors indeed found some evidence in 
favor of the cone geometry. In this paper, by performing similar comparison, we provide 
additional, more quantitative support to the cone conjecture. In addition, we test the 
scaling ( |1.2| ) and observe the power-law, but no "wiggles" in any dimensions. 

It is convenient to reproduce here the relevant formulae from [Q, ^ adopting them to 
our coordinates i \2.3\j . The comparison of the string's spacetime with the cone is done in 
Euclidean signature. The transformation is achieved by defining Ie = it which takes the 
metric (^^) to 

ds^ = f{p{r))e^^dtl + e^^{dr^ + dz^) + p{rfe^^ dnl_^. (2.21) 

As usual in the Euclidean signature, the imaginary time must be identified with period 
equal to the inverse temperature, tE G [0,4:7r/{D — 4)e^~^]. 

The cone 

The Euclidean double-cone metric over the base 5^ x S^^^ reads 

dsLne = dR' + R^ {dx^ + sin2 xK^dr^) + ^l^dnl_.^ (2.22) 

where x is the angle of 5^ ranging within (0, vr); t = tE and therefore k = ^{D — 4)e^~^ 
is the surface gravity of the black string. 
Comparing ( 2.21| ) with ( p. 22 ) one gets 



i2(r-,z) = e^('^-^)p(r)V§3|, 



sinx(r,.) = VP-4)/(p(0)^A(...)-C(...)_ (2.23) 
Kp(r) 

Next, one rewrites the remaining components of the cone metric in the (r, z) coordinates, 
di?(r, zf + jy-^Rir, zfdxir, zf = cl dr"^ + c2 dz^ 2 c3 drdz, (2.24) 
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where 



cl = {drRf + 
c2 = {d,Rf + 



1 

D-2 
1 

D-2 



R\drX)\ 



c3 = {drR){d,R) + 



1 



:R\drx){d,x)- 



(2.25) 



L> - 2 

Comparing with (2.21) one concludes that the cone geometry will indeed approximate the 
regions of the black string's spacetime where cl = c2 = (^'^'^) and c3 = 0. 

As an additional test, one can compare the Kretschmann scalar curvatures of the 
numerical and the cone metrics. For the cone this scalar is given by 

1 



K = RabcdR''^'"' OC 



(2.26) 



where the proportionality constant depends on the dimension (in 5D it is 48, in 6D it is 
72, in 7D it is 320/3 etc). For the metric (|2^ ) the expression for K is cumbersome but it 
is straightforward to obtain. 



The scaling and the wiggles 

In order to test the scaling (|1.2| ) one examines the behavior of a characteristic length 
scale of the nearly pinching strings as a function of some parameter along the NUBS 
branch. The characteristic variable whose behavior we test is Rmin- Identifying a good 
parametrization^ appears to be a difficult task in our case, since many potential parameters, 
such as k or S, for instance, are too noisy at the near merger limit due to numerical 
errors. In this paper, we use two different parameterizations: (i) By pi = Lhory which 
monotonically grows and seems to tend to a finite value, L^^^ , at large nonuniformity; 
(ii) By p2 ^ -8(7^, which is our low-level parameter, that apparently tends to zero at the 
merger. While the first parametrization is continuous across the merger point^, it is not 
obvious that it is also analytic at this point. As for the second parametrization - it becomes 
not analytic at the merger. The scaling ( |1.2| ) will be confirmed if the graph of log(i?772m) 
against \og{5pi) for i = 1 and i = 2 will be on average linear with periodic wiggles about 
it. The slope of the linear fit to the curve will provide 7i and the period of the wiggles 
about the fit is directly measured. 

3. Numerical implementation 

We fix the asymptotic length of the compact circle to the critical value L = = 2t: /kc, 
where kc is the critical Gregory-Laflamme wavenumber (see table 1 in [^), and generate 
the NUBS solutions by varying the constant Bq = B\^^i/2 ( p.lO| ). 

*A good parameter should be defined on both sides of the merger, i.e. along BHs and NUBSs branches, 
and it should be analytic at the merger point. In this case, the scaling exponent would not be affected by a 
transformation to another parametrization, provided the transformation is analytic. In the critical collapse 
context this is related with the notion of "universality", i.e. independence of parametrization. 

^On the BHs side we can define "Lhor" as a (normalized) proper length of the compact circle along the 
symmetry axis r = and between the poles of black hole along the horizon. 
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The actual solution of the elliptic equations (^^) is obtained using relaxation. To this 
end the equations are discretized on a lattice covering the domain of integration {{r,z) : 
< z < L/2,0 < r < Va}, where is the location of the outer boundary. The equations 
are written using a finite difference approximation (FDA) which is second order in the grid 
spacings. The relaxation of the non-linear equations ( |2.4D is achieved by incorporating 
a Newton iteration. Specifically, the fields at an interior grid point are updated 

according to 



old 



(3.1) 



d£^/d^{i,j) 

for ^ = A,B and C; £\fi{i,j) is the FDA equation of motion for ^ and c<j is a numerical 
factor. In the basic Gauss-Seidel relaxation method uj is set to unity. In our case we find 
that choosing uj = 1 for all fields causes divergence of the code. However, taking lob < 1, 
that is under-relaxing B, stabilized the scheme (usually we used ujb — 0.5). 

At the boundaries of the domain, boundary conditions (b.c.) are incorporated. While, 
for Neumann b.c. the update procedure is similar to ( |3.lD , for Direchlet or Robin (mixed) 
b.c. the equations are not solved. We use Direchlet condition ( 2.10| ) to update B at 
the horizon. The Robin boundary conditions are used at the outer boundary, where A 
and B are updated according to ( ^.11 ): A{z,ra) = A{z,ra — Ar)[ra/{ra — Ar)]^~^ and 
B{z,ra) = B{z,ra — Ar)[ra/{ra — Ar)]^~^, where Ar is the grid-spacing. C is updated 
according to 



(jr. 



(l-Mc)C°''(ra) + ^c 



C{z,ra-Ar) 



Ar 



old 



(3.2) 



where < fic < 1 is the "inertia" parameter. We generate most of our solutions using large 
inertia i.e. /Uc < 0.05. Using no inertia caused divergence. All other b.c. are Neumann. 

The relaxation begins with specifying an initial guess for the fields. Then, we sweep 
the grid using the red-black ordering and update the fields according to ( p.l| ) . The sweeps 
are repeated until the residuals drop below certain predefined tolerance (typically, 10~^). 

The hierarchy of meshes 

One of the disadvantages of relaxation methods is their relatively slow convergence. A 
well-known technique intended to accelerate the convergence is the full adaptive storage 
multi-grid method, see e.g. [|3l[. Unfortunately, in our case we could not find convergent 
multi-grid implementation^^ . 

The method that we found to work well makes use of several grids but in a different 
from the full multi-grid fashion. We begin with some initial guess on the coarsest grid and 
completely relax the equations using the algorithm described in the previous subsection. 
The solution is then interpolated into a finer grid where it is regarded as an initial guess. 
The equations are relaxed on this grid too, and the solution is passed to the next finer grid. 
The procedure is repeated on a desired amount of grids. The one-way manner that we 



^"it doesn't blow up, but the residuals just stop decreasing below certain limit. See appendix ^ for 
further discussion of the phenomenon. 
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Figure 2: The grid liierarchy, 3 mesiies are sliown. While the coarsest mesh covers the entire 
domain, the finer grids extend only to half size of the previous coarser mesh. The displayed field B 
(as well as A and C) is most variable near the horizon, justifying this "AMR" -type construction. 



propagate the solution along the grid hierarchy is to be contrasted with the bi-directional 
(V-cycle) communication between grids that is utilized by full multi-grid methods. 

We construct the mesh hierarchy using 2 : 1 refinement. Furthermore, in this sequence 
the radial extension of a finer mesh is taken to be half the extension of the previous mesh. 
The truncation is basically done because the high resolution delivered by dense grids is 
only necessary near the horizon where the fields are most variable^^. Note that number of 
the grid points in radial direction is conserved in this hierarchy. The boundary conditions 
used on finer grids do not change, except at the outer boundary, where we keep the fields at 
the values obtained on the coarsest grid. We verify that no mismatch between our NUBS 
solutions on any two subsequent grids arises along this boundary. Figure ^ illustrates the 
hierarchy and the field B on it in 6D. 

The initial guess on the coarsest mesh requires some caution. A priori we don't have 
a natural choice, so we use the simplest A = B = C = substitution. We found that this 
guess is relaxed effectively only for Bq below certain value (which roughly corresponds to 
mildly deformed strings) for larger i?o's it leads to a very slow convergence or to diver- 
gence. However, using a solution obtained at smaller Bq as an initial guess for higher -Bo's 
prevented the code from crushing and improved the convergence rate. 

Despite using this trick, on each particular initial grid we were able to find solutions 
only up to certain maximal Bq. For larger Bq^s the method diverges no matter how we 
tune up the parameters. In this case, we could continue to larger Bq^s by increasing the 

Hence our refinement is in tune with the adaptive mesh refinement (AMR) approach. 
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density of the initial mesh. A significant gain in the convergence rate was achieved by 
feeding the available lower resolution solution as an initial guess. 

As a result, after finding the NUBS branch we end up with a family of several initial 
meshes having different resolutions, and the multi-grid hierarchy attached to each one 
of them. We verify that the physical measurables obtained on various meshes within 
the hierarchy converge quadratically. When we had a family of several initial grids with 
overlapping parameters we verified the convergence rate as a function of the resolution on 
the initial mesh. Figures and |5| are made in six dimensions where we have such an 
overlap. The convergence rate of the quantities shown in figures is nearly second order. 

A verification of the constraints as well as further technical details of the numerical 
procedure are found in appendix ^ 

Overall, we find that the method performs very efficiently in lower dimensions and small 
nonuniformities. For large nonuniformities the convergence slows down and, in general, 
violation of the constraints grows. We could overcome the problems and accurately get the 
solutions with stronger horizon deformations by increasing the resolution of the numerical 
lattice. We had to stop when the mesh-size made the relaxation time unreasonably long. 
In higher dimensions the problems are more pronounced and begin earlier. In addition, 
the range of the parameters ub and fic, for which the relaxation converges, shrinks. In 
practice, using our current method we were only able to get the NUBS solutions in the 
dimensions between six to eleven. The following section summarizes our findings. 



4. Results 



We construct the branch of non- uniform black strings in dimensions 6 < D < 11 and 
extend it deep into non-linear regime. In 6D we find good agreement with previous results 
p!^ , in the range of A's obtained in these works, A ~ 4 in and A ~ 6 in and 
continue the branch up to A ~ 8. Below we describe the behavior of the characteristic 
physical variables and their dependence on the dimension. 

The typical behavior of the (normalized) entropy ( 2.16| ) in the verified range of the 
dimensions is depicted in figure ^. The key feature here is that the relative entropy grows 
with A, so that NUBSs are more entropic than the uniform critical string. In D < 9 
we find that at large A's the entropy asymptote to certain limiting value. These limiting 
values are listed in table |l| and they estimate the entropy of the merger string. In higher 
dimensions our numerical method does not allow extending the branch that far, in this 
case the numbers in the table are only lower bounds on the asymptotic /S"s. We find that 
the S{X) curve steepens as D grows and the saturations occurs earlier. For example, while 
in 6D this happens around A ~ 3, in 9D S reaches saturation just above A ~ 1.5. The 
entropy of the nearly pinching strings decreases with the dimension. 

We note that the S curve has the back-bending reported in |19|. However, the variation 
of S in the bending, 6S/S, is less than a percent in our range of A's. In addition, from 
the figure we learn that the specific value of A^;, beyond which the curve bends depends on 
the resolution of the numerical lattice used to find the solution, in a way that an increase 
in resolution increases Afef,. In 6D (shown in figure) we use a sequence of three resolutions 
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Figure 3: The normalized entropy (2.16) of NUBSs in 6D as a function of A. Shown are the 
solutions obtained using three mesh resolutions: low (pluses), medium (crosses) and high (stars), 
obtained by halving the grid-spacings. In the region where all three overlap we observe nearly 
second order convergence. The NUBSs are more entropic than that of the critical uniform string. 
In the limit of large A the entropy reaches saturation. The bending that shows up beyond certain 
A at lower resolution, at a higher resolution moves to larger A's in the manner that doesn't form 
clearly convergent sequence. 



with grid spacings /i, h/2 and h/A. The resulting X^bS are approximately 3.5,4.5 and 5.5, 
which doesn't seem to form a convergent sequence. Similar behavior is gotten in higher 
D's as well. 

The typical dependence of the temperature ( 2.15| ) on A is depicted in figure ^ which 
shows that the NUBSs are cooler than the critical GL string. The temperature decreases 
with A and in 6 < -D < 9 it is seen to tend to certain limiting values. In higher dimen- 
sions our numerics loses stability before the saturation is achieved. The minimal attained 
temperatures are listed in table |l|. The features such as the steepening of the T{\) curve 
in higher dimensions and the bending are similar to those observed for the entropy. The 
temperature of the nearly pinching strings in various dimensions increases with D. 

In figure ^ we plot the normalized mass as a function of A. The mass is computed in 
two ways: (i) From the asymptotic charges according to ( 2.13D , the result is designated 
by markers; (ii) Integrating the first law ( |2.18| ) according to (|B.3| ), this is represented by 
solid line. The figure demonstrates that both methods yield comparable results and the 
agreement further improves as we increase the resolution. In all verified dimensions the 
mass along NUBS branch is higher than that of the critical string. All other features of 
the mass curve are similar to those of the entropy curve. At large A, in 6 < D < 9, the 
mass reaches saturation, the asymptotic values recorded in table ||. In higher dimensions 
the table shows only the maximal masses obtained rather than asymptotics. The mass of 
the nearly pinched strings is a decreasing function of the dimension. 
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Figure 4: The temperature in 6D as a function of A. Same conventions as in figure ||. 

6D 




Figure 5: The mass as a function of A in 6D, same conventions as in figure 0. The markers 



designate the mass computed from (2.13), the solid hues indicate the resuh obtained by integration 



of the first law, see (B.3). Both methods agree well and the agreement improves at higher resolution. 



Smarr's formula ( 2.17] ) relates the local horizon variables S and T with the asymptotic 
ones, m and r, and as such it is an important indicator of global accuracy of a numerical 
method^^. In figure ^ we plot the ratio between the two sides of ( |2.17| ) in six and eight 
dimensions. We observe that in general the formula is well satisfied. The scattering of 



^■^It turns out that the tension r along the NUBS branch cannot be extracted accurately enough in our 
numerical implementation. However, this fact does not have major impact on the accuracy of the Smarr 
relation. See appendix Ira for analysis and discussion. 
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Figure 6: The Smarr formula ( p. 17 ). Shown are the ratios between the two sides of (2.17) in 6D 
(squares) and 8D (stars). In each dimension, there are three sequences of points obtained at three 
numerical resolutions, with the sequence at largest A corresponding to the largest resolution. In 
either case, the points concentrate near 1 such that the violation of Smarr's relation never exceeds 
the 10% level. 



the points does not exceed 10%, a large fraction of which is probably due to inaccurate 
r. For a given resolution the accuracy decreases for progressively non-uniform solutions. 
However, the accuracy improves at higher resolutions, indicating convergence. 

To analyze the geometry of the NUBSs we plot in figure the minimal and the maximal 
areal radii of the horizon and its proper length along the compact circle as a function of 
A. We find that while Rmin shrinks monotonically with A, Rmax and Lhor asymptote to 
finite values in the large A limit. The behavior of the curve Rmaxi^) is again similar to the 
behavior of the entropy. The LhorW curve is similar too, with the only difference that it 
does not have any back-bending. As the dimension grows the asymptotic values of Rmax 
and Lhor are getting smaller, see table |^. 

Figure |8| shows the embedding ( |2.19| ) of the black string horizon into fiat space. The 
maximal and the minimal radii of the horizon are explicitly seen to occur at z = 0,L and 
at z = -L/2 respectively. The proper length of the circle along r = const is maximal at 
the horizon but for large r, L ^ 1 exponentially fast. Note that the transition from a 
uniform to a non-uniform black string causes an expansion of the compact circle near the 
horizon. In this sense we have an Archimedes effect for nonuniform black strings. The 
effect weakens in higher dimensions, see table |l[ 

The local cone 

We compare the local geometry in vicinity of the "waist" of strongly nonuniform strings 
with the cone metric. The relevant spacetime region is illustrated in figure |. The cone 
sketched in the right panel is predicted to approximate the local merger geometry of the 
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Figure 7: The typical beiiavior of tlie geometric variables: i?mm decreases steadily with A, R,, 
and Lhor asymptote to finite values. Same conventions as in figure 0. 



D 




m 


T 


S 


^max 




Lhor 


6 


8 


2.26 


0.73 


3.3 


1.75 


0.103 


1.45 


7 


4.45 


2.37 


0.793 


3.25 


1.57 


0.16 


1.43 


8 


4.1 


2.28 


0.86 


2.83 


1.42 


0.156 


1.34 


9 


1.7 


1.98 


0.91 


2.3 


1.32 


0.3 


1.23 


10 


0.46 


1.5 


0.96 


1.7 


1.22 


0.64 


1.07 


11 


0.2 


1.14 


0.99 


1.22 


1.13 


0.81 


1.04 



Table 1: The thermodynamic and the geometric variables at the maximal A to which we could 
extend the NUBS branch. In the dimensions between six to nine the behavior of the variables is 
similar to this plotted in figures |3[ ^, ^ and |^. Namely, in these dimensions the listed numbers are 
close to the asymptotics. {Rmm does not reach saturation but continuously decreases. The shown 
values correspond to the minimal attained Rmm-) In higher dimensions the shown numbers are not 
asymptotics, but only bounds on the true limiting values. 



marginally pinching string. In this case one would expect that a non-pinching, but strongly 
deformed string still carries signs of the cone geometry in the region bounded by the circles 
shown in figure. We will now demonstrate that with high precision the geometry near the 
"waist" is indeed cone-like. 

We explicitly compare the numeric and the cone metrics in six dimensions where we 
have our most non-uniform solutions and the highest resolutions. First, according to ( 2.23| ) 
we obtain R and x iii terms of (r, z). The extraction of x appears to be delicate |pT|| 
because x is obtained from sinx, which for our numerical solutions with finite A behaves 
somewhat differently from what it should for the cone. Specifically, at the axis z = 
that corresponds to x = "^72, the cone must have sinx = 1. However, since the cone apex 
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Figure 8: The embedding of tlie string's horizon into flat space ( 2.19 ). The right panel shows the 



cone and the region (between the circles) where the local geometry is compared with the cone. 



is missing at finite A, sinx obtained from (2.23) happens to be greater than 1 at small r's, 



see figures |^ and |^. At the same time, figure ^ shows that along the axis sine is close to 1. 
The peak in the curve, that appears because the cone's tip is missing, drifts to the origin 
as A grows indicating that the cone apex becomes resolved better in the large A limit. 
Clearly, having sin x\axis > 1 obstructs the extraction of x- A practical remedy to the 



problem was suggested in |21], where the sine was "renormalized" by its value at the axis, 



forcing it to the correct, sin Xaxis = value. The renormalization is achieved by 

■ ^ M sinx{r,z) 

Sm x{r,Z)\ corrected = —. , (4.1) 

sill Xaxis 

where in order to normalize outside of the axis one finds sinx at the axis and proceeds 
along R = const curves using this value. 

Having obtained R and we determine the metric components ( p.25 ) cl,c2 and c3 



and plot them in figures [10| and |11| for A ~ 6.3 and A ~ 8 respectively. There is an excellent 
agreement with the cone prediction. It is required that cl = c2 = e'^ ^ and from the figures 
we learn that while c2 and cl differ from e^^ by several tens of percents far away from 
the waist, the match grows to a level of more then 90% closer to the corner. Expectedly, 
the agreement worsens again in the immediate vicinity of the apex, see figure H As for c3, 
which must vanish for the cone - it is two orders of magnitude smaller than c2, cl in the 
relevant region, and this is indeed consistent with the prediction. 

Another test involves comparison of the Kretschmann curvature invariants computed 



for the numerical and the cone metrics, see figure 12. The plot indicates that the invariants 



are comparable around the "waist", where the mismatch between them is only about 10%. 
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Figure 9: sin^ along the z = L/2 axis as determined from (2.23) in 6D. For the cone sinx|a£cis = 1, 
but for a finite A NUBS this is not so. Since the cone tip is missing, the curve initially climbs steeply 
to values greater than 1, stays close to unity and finally deviates from one farther away from the 
waist. For increasingly nonuniform strings the peak tends to the origin and the entire curve gets 
consistently closer to unity. This indicates that the cone tip is closer approached. The markers in 
one of the graphs designate the actual resolution of the numerical lattice. 



Next we verify the scaling ( |1.2| ). A characteristic quantity with dimensions of length 
whose scaling we will test is the minimal areal radius of the horizon Rmin- It is measured 
directly from the metric and hence it is the cleanest, from the numerical point of view, 
variable. One of the parameterizations of the NUBS branch that we use in this paper is 
given hy p = Lhor- The horizon length apparently remains finite at the pinch-off, and as 
we saw it doesn't exhibit any bending. Assuming that Lhor is growing monotonically we 
obtain the limiting length, p* = L^^^ by extrapolation. 

Figure 13 plots log(i?mm) as a function of 5p = log(LJ^°^ — L/ior) in six, seven, eight 
and nine dimensions. Figure shows that for small p's the scaling of Rmin is not a power-law. 
However, the behavior is approaching a power-law dependence in the p ^ p* limit. This 
is most pronounced in 6D where we have the longest data series. The "runaways" of the 
data points that appear at the end of each sequence with a given resolution is probably 
related to the loss of accuracy. 

The insets show several points that correspond to the largest non-uniformity. These 
points can be fit quite accurately by a straight line (shown), and this indicates the power- 
law scaling. The corresponding exponents are given by the slope of the fits and they are 
summarized in table ^ These are only approximate values because of the numerical errors 
related to the fitting procedure (whose accuracy is within 5%), and other manipulation 
such as extrapolation that we used to find p* etc. Besides, the values in table should 
strictly speaking be regarded as upper bounds on the exponents, rather than the exponents 
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Figure 10: 6D; A ~ 6.3. The metric components ( 2.25 ) tliat must satisfy cl = c2 = and 



c3 = if the cone indeed approximates local geometry near the "waist" which is located at the 
bottom right corner, cl and c2 agree very well with the cone prediction: we see approximately 
90% match slightly away from the corner. c3 is two orders of magnitude smaller than c2 and cl, 
that is consistent with the cone prediction as well. The sub-domain shown here is covered with the 
numerical lattice of approximately 200 x 200 mesh-points. 



themselves. It is hard to say how close to the actual 7's the bounds are, but at least in six 
and seven dimensions they seem to provide a good estimate, see figure 

In principle, we could parameterize the NUBSs branch differently, for instance, by k or 
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Figure 11: 6D. A ~ 8 and we use same conventions as in figure but remove the apex. The 
agreement with the cone prediction is somewhat better than in figure nOl 



by S. In practice however, all measurables except i/ior have the back-bending that hampers 
determination of their asymptotic values, which are crucial for forming 5p. On the other 
hand, the control parameter, Bq^ , that we use to generate our solutions is a monotonic 
function of A, and it naturally parameterizes the NUBSs branch. A disadvantage of this 



parametrization is that it is not defined on the BHs side. Figure 14 depicts log(i?. 



vs logi?Q . The behavior is not a power-low for mildly deformed strings, however as 
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Figure 12: 6D. A comparison of tlie Kretsclimann invariants computed for the cone and the 
numerical black-string with A ~ 6.3. A good match is in the vicinity of the "waist", which is 
located at the right bottom corner of the figure. 



D 


6 


7 


8 


9 


7 in Lhor 


0.6 


0.54 


0.41 


0.33 


7 in B^^ 


2.1 


1.8 


2.2 


1.4 



Table 2: The scaling exponents 7 in Rmin oc Sp^ in various dimensions for two parameterizations. 
The first row shows the upper bounds on 7 in the parametrization by Lhor- The second row is the 
lower bound on 7 in the parametrization by Bq^ . Distinct 7's for a given D reflect the fact that 
the transformation between the parameterizations is not analytic. 



indicated by the straight lines shown in figure, the dependence becomes such in the near- 
merger regime. The slopes of the lines provide lower bounds on the critical exponents 7, 
we list these in table § and observe that 7 ~ 2 in all recorded dimensions. Comparing 7's 
in both parametrization we see that they are distinct, therefore the parameterizations are 
not analytically related. 

Having estimated the exponents we can ask about the fine structure, or the "wiggles", 
that according to the prediction should appear as a wavy structure about the linear fit in 



figures |1^ and |lj. By examining figure |13| one infers that in six and seven dimensions there 
are no any wiggles in the range of the achieved nonuniformity, and in higher dimensions 
the accuracy of the data is insufficient for deciding about existance/absence of the wiggles. 



From figure 14 it follows that the wiggles do not appear in any of the plots. 



5. Discussion 



Our numerical construction of the branch of the NUBS solutions confirms the trends ob- 
served in the perturbative analysis of and lends further support to the phase diagram 



- 21 - 



6D 



-0.5 



-1.5 



-2.5 



o low 
o medium 
o high 
linear fit 






-0.2 




-0.4 




-0.6 




-0.8 


c 
E 


-1 


CC 




-1.2 




-1.4 




-1.6 




-1.8 




-2 



o low 
o medii 
o high 



o o 
o o 



-3.5 



oo 
o o 



-3.4 -3.1 -2.1 

-3 -2.5 -2 -1.5 -1 

log(L™'<-L ) 
hor hor' 





y = 0.33*x + 2.5e-05 



o low 
high 



Figure 13: The scaling of log(i?mm) as a function of Sp = log(I/™J^^ — Lhor) in various dimensions. 
We show the sequences of data points obtained at several resolutions, with the largest resolution 
sequence corresponding to the smallest Sp. The inset shows only the highest resolution data that 
corresponds to the most deformed strings. The solid line, which is the linear fit, is seen to approxi- 
mate the data points well. This indicates that at large A's we have a power-low. The corresponding 
exponents are given by the pre-factors in the displayed fitting equations. 



|T]. We managed to extend the solutions branch into strongly nonuniform regime and this 
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Figure 14: The scaling of log(i?mm) as a function of Sp = \og{BQ ) in various dimensions. Same 
conventions as in figure In the limit of large deformations the scaling is clearly a power-law as 
indicated by the shown linear fit. 



allowed us to estimate the physical parameters of the nearly pinching strings, see table |I|. 
We expect the estimates to be good in 6 < D < 9, where we could approach the merger 
closely, and they are only bounds (upper or lower, depending on the variable) in higher 
dimensions. 

It follows that when the spacetime dimension grows the mass and the entropy of the 
near-merger solution decrease and their temperature increases. This is in accordance with 
the expectation Q that in higher dimensions the trends in the behavior of the (normal- 
ized) thermodynamic variables are reversed, see diagram ||. This is best illustrated by the 
temperature, which is expected to change its trend in D > 12. One notes that while in 
6D the limiting temperature along the branch is smaller than 1 by about 30%, in IID it is 
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smaller than 1 by less than 1.5%. (Of coarse, in this dimension we have not extended the 
branch into really asymptotic regime. However, the shape of the curve T(A) suggests that 
the asymptotic value is not much smaller than the estimate.) 

It is interesting to note that our numerical method performs well for weakly deformed 
strings in D < 11, but in 12D (and higher) it fails to converge for any, no matter how 
small, deformations. Recalling that the trends in the temperature are expected to change 
at about D ~ 12 we tend to attribute this behavior to the fact that due to limited accuracy 
our numerics is incapable of deciding on the sign of T — 1, even if it can determine correctly 
other thermodynamic variables. This probably leads to a conflict that eventually results 
in the divergence. 

Examining our most deformed strings we provide strong quantitative evidence that the 
local geometry in vicinity of their waist is cone-like. The deviation from the cone geometry 
is less than ten percent, which leaves little room for any solution other than the cone. The 
cone's presence is further (indirectly) supported by appearance of the power-law scaling 

(HI)- 

The power-law behavior (more precisely - its onset) of the characteristic length scale 
Rmin is observed and depicted in figures |l^ and ^ The estimates of the corresponding 
critical exponents are collected in table ||. We used two parameterizations that are not ana- 
lytically related, and hence the corresponding 7's are not equal. While the parametrization 
by Bq^ is not defined beyond the merger point (along the BHs branch), the parametriza- 
tion by L/joj. can be made continuous across this point (though it is not obvious if it is 
analytic at the pinch-off); in this sense Lhor parametrization is superior. 

Proving that L/ior is analytic at the pinch-off, will qualify the parametrization as "good" 
(see footnote 8). In this case^^ we can compare the estimated 7's with the prediction of |po|| , 
where it is proposed that jth = 1/4 in any D < 10. Clearly, our estimates differ from 1/4 
and there can be several reasons for this. First, in general, one should not really expect for 
the actual exponents to coincide with those obtained in the linear analysis of [^0|. Second, 
our estimates may be too crude and the more accurate 7's will indeed be 1/4. While this 
is possible, we believe that at least in six and in seven dimensions our estimates cannot be 



significantly different from the actual exponents, see figure 13. Finally, it is likely that the 
truth lies somewhere in the middle, and so are the actual values of the critical exponents. 

In conclusion, we note that there are still many aspects of the phase diagram || awaiting 
for a confirmation. This includes a construction of the black-hole branch in higher than six 
dimensions and of the black-string branch in D > 11. Despite the progress in understanding 
the pinch-off geometry reported in this paper, one still might wish to approach the merger 
closer and to do so also form the "black-hole side" . Our current relaxation method seems 
inadequate for achieving these goals due to its slow convergence and instability at the 
strongest non-uniformity. We believe that in order to address the mentioned issues our 
scheme will require major revisions (e.g. finding a stable multi-grid implementation), or 
one would probably have to come up with a completely different approach. 



^If Lhor won't happen to be analytic at the merger, the exponents are incomparable. 
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A. Uniform black string in conformal coordinates 

The metric of the uniform black string in Schwarzschild coordinates reads 



ds' = dsi.,^^^ + dz' = -f{p)dt' + fipr'dp' + dz' + p'dni^^ 
f{p) = l-{p,/pf-\ 



(A.l) 



where pQ designates the horizon location. 

In order to bring the (r, z) part of the metric into the conformal form we change the 
coordinates such that 



ds'^ = -f{p{r))dt^ + dr^ + dz^ + p{rfdnl_^ 
dr/dp = r^'\p) 

r[p) = ro + ' 



(A.2) 



where tq is the integration constant, that designates the location of the horizon in the 
conformal gauge (without loss of generality we choose ro = 0.) The integration gives 



1 



-Po 



d-5 ^ -1-/^2-^1 [ (i_3i 2' d-3' 



r( 



2{d-S) 



d-3 



for d = A; 



for d 



, for d > 5, 



(A.3) 



where is the hypergeometric function. 

Our numerical implementation in the conformal coordinates requires knowing p{r). 
Explicit inversion of ( [A.3| ) exists only in d = 5, in any other dimension the inversion must 
be achieved numerically. 



B. Details and tests of the numerics 

In this appendix we provide further details of the numerical procedure, verify the con- 
straints and describe how the asymptotic charges are obtained. 

Table ^ summarizes the measures of the domain of integration and the sizes of the 
coarsest meshes in various dimensions. In order to minimize the mesh size (and, as a 
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iterations/10 



Figure 15: A typical behavior of the residuals for several sets of the parameters lub and ^c- In 
some cases the residuals blow up after finite amount of iterations or show very slow convergence. 
In other cases the convergence rate is faster. 

result, the relaxation time) without compromising accuracy we choose the location of the 
asymptotic boundary ra as close to the origin as possible. To achieve this, we try several 
r^'s and select the minimal one for which the thermodynamic variables S and T undergo 
less than 0.5% variation if is increased. 



D 


L/2 


ra 


low 


medium 


high 


6 


2.4758 


7 


201 X 65 


401 X 129 


801 X 257 


7 


1.9884 


6 


121 X 41 


241 X 81 


481 X 161 


8 


1.6982 


5 


101 X 35 


201 X 69 


401 X 137 


9 


1.5032 


4.5 


91 X 33 


181 X 65 


361 X 129 


10 


1.3659 


4 


91 X 33 


181 X 65 


361 X 129 


11 


1.2566 


3.75 


91 X 33 


181 X 65 





Table 3: Location of the outer boundary ra and sizes of the coarsest meshes used in various 
dimensions. L/2 = n/kc, where kc is the critical GL wavenumber, see table 1 in 

Relaxation parameters 

The parameters lob and fic are selected after several trial runs with similar geometry and 
resolution but with different parameters. Figure ^ shows the behavior of the L2-norm of 
the residuals in several cases. (The L2-norm of matrix M on mesh [1 . . . Nr] x [1 . . . Nz] 
is defined by ||M||l2 = ^JYlij \^ij\/^rNz.) While for some sets of the parameters, the 
residuals blow up or converge very slowly, for other sets the residuals decrease faster. 
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We have not attempted to locate the "optimal" parameters that might accelerate the 
convergence rate. 

The necessity to use under-relaxation for B, strong asymptotic inertia for C and to 
fine-tune the corresponding parameters into a narrow range is bewildering. In addition, the 
convergence/divergence of the code is strongly affected by the initial guess. As described in 
section ^, for the relaxation to converge one has to choose the initial guess "close enough" 
to the final solution, otherwise the code diverges. 

We tend to attribute this unhappy situation to the fact that (1 — e'^^~'^^)/p'^ terms 
in (|2.4| ) have "wrong sign" (see also the discussion in Q). In order to explain what is 
wrong with the sign, we linearize around some solution. Doing so for ^' = i? or C we 
end up with the equation of the form V^^* + m^"^ + . . . = 0, where m depends on the 
background solution and the position, and the ellipsis designate terms containing first 
derivatives. Ignoring for a moment these terms we locally get a Helmholtz-type equation. 
A boundary-value problem for this type of equation is ill-posed if m is smaller than the 
mesh box, namely if the tachyon defined by the equation fits into the box and might get 
exited. This clearly destroys any hopes to relax such an equation numerically. 

As for the original set of equations (|2.4| ), the finite range of the fine-tuned parameters 
ub and /xc and the sensitivity to the initial guess apparently indicates that while outside of 
this range the "wrong-sign" causes divergence, inside this range the tachyon is suppressed, 
and the relaxation converges. We note also that for increasingly non-uniform strings and 
in higher dimensions the convergence window shrinks (in lOD, for example, fic < 0.01) 
The sensitivity to specific values of ujb^^^c and to the initial guess may explain the failure 
of the full multi-grid technique. The method operates on grids of variable size and this 
probably causes the excitation of the tachyon, since it may happen to fit into one of the 
grid boxes. 

Note that at least one of the wrong-sign equations is a direct consequence of working 
in a higher dimensional spacetime, D > 4, that has 0(A^-|-1) isometry with D — 3>N>2. 
Indeed, in this case the most general metric is given by 

ds'^ = dsl_j^{X) + e^^^^Un%. (B.l) 

Dimensionally reducing the Einstein-Hilbert action over the sphere and varying it to get 
the equations (see for example appendix A of Q where this is elaborated for N = D — 3), 
one finds out that C is governed by 

AC-G(5"'5)e-2^ + --- = 0, (B.2) 

where ga/3 is a metric on X, A = gx^'^'^a{gx^g°^^dp) and ellipsis designate terms with 
first order derivatives and second order derivatives of gaj3- The function G{g°^^) depends 
on the metric 5"^, and it is positive in D > 5 and zero otherwise. The equation (|B.2| ) 
has wrong sign in the sense described above, and this property is invariant under gauge 
transformations of the manifold X. 
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Figure 16: V and rel V for a strongly nonuniform string (A ~ 3.5) in 6D. Tlie maximal deviation 
from zero occurs at the horizon, near z — L/2. Outside this location the constraint is small. The 
relative constraint is small too. Altogether we conclude that the constraint is well satisfied. 



Constraints 

A violation of constraints is inevitable in any numerical solution. Here, we examine the 
convergence rate of the constraints U and V as a function of the grid-spacings. In addition, 
we verify that the (properly normalized) constraints are small. 

Usually we solved the equations on the three-grid hierarchy. We computed L2-iiorms of 
both constraints on these grids and checked that typically the ratio (||Z//||'' — ||Z^||'^/^)/(||Z//||''/^- 
||Z//||'^/^) > 3 and the similar ratio for V is slightly larger than one. This indicates that while 
the constraint U converges well, the convergence of V is poor. We tend to attribute this 
behavior of V to a growth of the discretization errors near the horizon due to the terms 
containing the singular l/vT- We will demonstrate now that V is actually very small 
everywhere except in the horizon vicinity. 

In order to get an idea of how well the constraints are satisfied, we define the "relative 
constraints": rel U = U/^\teims mU\ and rel V = V/^|terms inV|, where the sum 



is over the absolute values of all additive terms forming lA or V. In figure 16 we depict 
V and its relative value. V is maximal near the horizon around z = L/2 and it is small 
everywhere else. On the other hand the relative constraint is small at this corner, but it 
is large along the remaining section of the horizon (and its neighborhood), where the V 
itself is small (typically of order of 10~'^ — 10~^). Therefore, both plots are complimentary 
and together they indicate that the constraint is well satisfied. The peak that V has at the 
horizon, see figure, gets sharper at higher resolution, but at the same time it also moves 
towards the corner. This is why the L2-norm of V showed poor convergence. 



The second constraint lA and its normalized counterpart are illustrated in figure 17. lA 
is small everywhere and it consistently diminishes with increasing grid density. The relative 
lA is very small as well, except near the outer boundary, where it gets large simply because 
lA is extremely small there and thus the ratio that forms rel lA becomes of order unity. In 
summary, the plots are complimentary and indicate that the constraint is well satisfied. 
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Figure 17: The oc U and rel U constraints for a strongly non-uniform string (A ~ 3.5) in 
6D. The constraint itself is essentially small, with the maximum reached near the horizon. The 
relative constraint is also very small. The observed growing towards the outer boundary is caused 
by smallness of lA, see the insertions. Both plots are complementary and together they indicate 
that the constraint is well satisfied. 



Extraction of the charges 

The mass and the tension can be obtained from the asymptotic constants (2.11) using 



(|2.12 ). The constants are extracted by fitting the fields near the outer boundary with 



analytic expressions of the form ( |2.11 ). It turns out that doing the fitting near is 



unreliable. Hence, in practice we compute the (averaged along Tq) fields and their r- 
derivatives, and use this to integrate the fields from ra to a larger (typically ^ 10 ra) 
distance. The constants are then extracted from fitting near this boundary. We integrate 
( |2.4| ) to get A and C. In order to find B we solve the first order equation obtained after 
the second derivative d^B in ( p.4[ ) is eliminated with the aid of the constraint Gzz- (Using 
the second order equation ( |2.4| ) causes blow-up, see also |l^). Of course in the asymptotic 
equations we drop the (negligible at large r's) terms containing z-derivatives. 

We find that the constant a is essentially robust as it has very little dependence on the 
details of the integrating-fitting procedure. On the other hand, b is very sensitive to these 



details and might vary as strong as by 50% when some of them change (see [l^i IMl 
related discussions). Therefore the mass, which is essentially determined by a, is obtained 
accurately, but the tension, which is mainly contributed by 6, is not, see ( 2.12D . In fact we 



could only confirm that the relative tension is consistently within the range [0, \/{D — 3)] 
in dimensions D = 6, 7, 8. 

On the other hand, the local horizon quantities S and T are more accurate since they 
are measured directly without any additional manipulations that can enhance numerical 
errors. Hence, it makes sense to try computing the mass by integrating the first law ( 2.1^ ) 
along the solutions branch, 

mi{\) = m(0) + r(A') — d\'. (B.3) 
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We found that the mass computed in this way is close to that obtained from the asymptotic 
constants. The agreement improves when the numerical resolution increases, see figure 
|5|; this indicates the convergence of the numerics. The discrepancy between the mass 
computed in either way is never above 5%. 

Encouraged by this and by the apparent accuracy of Smarr's formula ( |2.17 ), shown in 



figure ^, one might hope to determine the tension without need to invoke the asymptotic 
integration-fitting procedure. Unfortunately this is not the case. From ( p. 17 ) it follows 
that the tension is given by the difference of two large numbers, 

T = {{D -3)M - {D -2)ST) /L (B.4) 

each of which becomes larger than tL by more than an order of magnitude for A > 0.1. 
One has to know M, S and T accurately enough to be able to reliably extract f from ( |B.^ ). 
It follows that in our case f computed from (B.4) is very noisy, sometimes with more than 



50% scattering. In addition, the accuracy worsens in higher dimensions, partly due to the 



D-dependent factors in (B.4). Consequently, this method does not produce more accurate 
estimate of r. 
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